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Abstract. We propose an extension of the semi-analytical solutions derived by Lin et al. (1965) describing the two- 
dimensional homologous collapse of a self-gravitating rotating cloud having uniform density and spheroidal shape, 
which includes magnetic field (with important restrictions) and thermal pressure. Fhe evolution of the cloud is 
reduced to three time dependent ordinary equations allowing to conduct a quick and preliminary investigation 
of the cloud dynamics during the precollapse phase, for a wide range of parameters. We apply our model to the 
collapse of a rotating and magnetized oblate and prolate isothermal core. Hydrodynamical numerical simulations 
are performed and comparison with the semi-analytical solutions is discussed. Under the assumption that all cores 
are similar, an apparent cloud axis ratio distribution is calculated from the sequence of successive evolutionary 
states for each of a large set of initial conditions. Fhe comparison with the observational distribution of the starless 
dense cores belonging to the catalog of Jijina et al. (1999) shows a good agreement for the rotating and initially 
£f>l ' prolate cores (aspect ratio ~ 0.5) permeated by an helical magnetic field (~ 17 — 20/iG for a density of ~ 10 

cm" 3 ). 
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^ ■ 1. Introduction and to give strong benchmarks to the numerical studies. 

However in spite of its importance, very few analytical 
The stars form through the gravitational collapse of dense tools are available to describe the dynamics f the precol- 



cores located in the molecular clouds. These condensations 



lapsing or starless phase. The present paper is an attempt 
are complex and involve gravity, hydrodynamics, magnetic ^ Q ^\ ^ ms g a p 



5^ 

■ and thermal processes 

, , , . ,. , , . In this paper, we first present new solutions of the 
Ine core dynamics has been investigated numerically 

, i ,i * iii i .1 -i i gravo-magnetic condensation. They are an extension of so- 

by several authors. Although the numerical approach un- ? ° . , . ^J.. . . . 

, , , ,. „ ,, , , , lutions obtained by Lin et al. (1965 for the gravitational 

doubtedly allows a deep understanding of the cloud dy- , . , . . , , , \ , , 

., , n , n i ;• i ; t n ; part and simultaneously by Aburihan et al. 2001 and by 

namics, it does not replace the analytical studies, first , , , TT , „ . , „ „ TT ^ ; TT 

, , , . , , . ,. ., , . ',. the author (Hennebelle 2001) for the MHD part. We then 

because analytical work gives more explicit descriptions \ . ' ^ 

, , , , , i r • i .1 i use these solutions to study the phenomenology of the col- 
and second because tne complexity ol numerical metnods 

, . . , . r , , . r lapse ol oblate and prolate, isothermal, magnetized and 

keeps increasing, making a full exploration of the space , , 

i,i , . i , ..J,, non magnetized clouds and compare the result with the 

parameters and tneir use by non numerical experts dim- . ° ..... „, 

evolution obtained with numerical simulations. The solu- 

\ , . ,, r i rr i i ! i- i tions allow to calculate the cloud shape evolution. With 

Analytically, most of the efforts have been dedicated .. . ln . n . , ., ,. r 

, , . .. n , . , , . . , the assumption that all cores are similar, distributions of 

to the description of the isothermal gravitational collapse . . ,. ,. . , T , r ,, 

„ TT Xi apparent axis ratio are predicted. We then compare these 

Larson 1969, Penston 1969a, 1969b, Hunter 1977, Shu ., .. .., v , . e 

„ T , . , „ „ X , , distributions with the apparent axis ratio distribution ol 

1977, Whitworth & Summers 1985, Bouquet et al. 1985). , , , X , ■ . , , c t - 

_, ' . , . . ..... . the starless dense cores belonging to the catalog of Jijina 

1 hese solutions which assume self-similarity are very use- ^ (1999) 

ful to get a physical hint on how the collapse proceeds 

Indeed the geometry of the dense cores is an important 

Send offprint requests to: P. Hennebelle, e-mail: issue since it depends on the physical processes that drive 

patrick. hennebelle@ens.fr their evolution. Previous investigations by Myers et al. 
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(1991) and Rydcn (1996) conclude that the observational 
data are more consistent with axisymmetric prolate 
cores than with axisymmetric oblate cores. The same 
conclusion is reached by the recent study of Curry (2002) 
with a different sample of cores and by Hartmann (2002) 
who shows that the major axis of the cores in the Taurus 
cloud is preferentially aligned with the filament in which 
they are embedded. Two recent studies, (Jones et al. 
2001 and Goodwin et al. 2002) analyse the observational 
apparent axis ratio distribution derived from the catalog 
of Jijina et al. (1999) and conclude that the core apparent 
axis ratio distribution is compatible with triaxal cores but 
being more nearly oblate than prolate. From a theoretical 
point of view, the equilibrium solutions of a rotating 
and magnetized cloud (Mouschovias 1976, Tomisaka 
et al. 1988) present an oblate shape except if there is 
a substantial toroidal magnetic component (Tomisaka 
1991) or if the structure of the poloidal magnetic field is 
such that it compresses the cloud (Stahler & Curry 2001). 
Curry (2000) also shows that two-dimensional solutions 
of the isothermal Lane-Emden equation present a prolate 
configurations. On the other hand, recent studies show 
(Nakamura et al. 1993, Fiege & Pudritz 2000a, 2000b) 
that prolate cores can be formed from a magnetized 
filament (with poloidal and toroidal field) through a com- 
bination of the gravitational and the magnetic sausage 
instabilities. Fiege & Pudritz (2000c) also demonstrate 
that the axis ratio of the prolate magnetized cores is 
compatible with the observational apparent axis ratio 
distribution. 

In the second Section of the paper, we derive the so- 
lutions and reduce the self-gravitating MHD equations 
to three time dependent ordinary equations. We discuss 
their mathematical and physical properties, their poten- 
tial and intrinsic weaknesses. In the third Section, we 
solve these ordinary equations numerically and describe 
the phenomenology of the collapse of an isothermal dense 
core. We also compare the semi- analytical solutions with 
the results of the hydrodynamical numerical simulations 
of an isothermal dense core. The fourth Section is devoted 
to the comparison between the theoretical and observa- 
tional apparent axis ratio distributions. A discussion and 
summary are given in the fifth part. 



through the dilatation groups (see e.g. Olver 1986 and 
Ferrara & Shchekinov 1996 for an introduction to the ap- 
plication of Lie groups to differential equations) . 

In this paper, on the contrary, we first assume a special 
and simple spatial dependence and put all the complexity 
in the time dependence of the fields. This technique has 
been successfully applied by Lin et al. (1965) to predict the 
eccentricity evolution of a cold cloud and by Hennebelle 
& Perault (2000) to study the thermo-magnetic conden- 
sation. 

2.1. The equations 

We consider the perfect MHD equations of self-gravitating 
gas with polytropic equation of state and assume an ax- 
isymmetric geometry. With the usual notations, we have 
in cylindrical coordinates: 



d t p + — d m {m Vvj p) + d z (V z p) = 0, 
w 



(1) 



P UtV„ + v„d m v„ + v z d z v m - = -dt.P 

+pd^ - —Bod^{mBe) + — B z {d z B m - 8^B Z ) , (2) 

fl W flQ 



VgV^ 



P d t V e + V„d a V e + V z 8 z V e + 



— B w d„{wBg) + — B z d z Bg, 

p w ho 



(3) 



P (dtV z + V w d^V z + V z d z V z ) = -0 Z P 

+pd z $ - —B d z Bg - —B m (d z B^ - d m B z ) , (4) 
Po Po 



dtiPp- 1 ) + V^{Pp-"<) + V z d z (Pp-^) = 0, 



1 



(5) 



(6) 



2. Homologous solutions of the gravo-magnetic 
condensation 

In this section we look for self-similar solutions of Eqs. |JI}- 
<|10|) and reduce this system into a system of three ordinary 
differential equations of the time. 

In most of the classic studies of self-similar solutions 
(Larson 1969, Penston 1969a, 1969b, Hunter 1977, Shu 
1977, Whitworth & Summers 1985, Bouquet et al. 1985), 
a special and simple time dependence of the fields is as- 
sumed leading to ordinary equations of the self-similar 
variable. These solutions are indeed related to the in- 
variance of the self-gravitating hydrodynamical equations 



-D^mB^) + 8 Z B Z = 0, (7) 
w 

d t B^ - d z {V^B z - V z B m ) = 0, (8) 

d t B e + d z (V z B e ~ V e B z ) - dMB m - V„B e ) - 0, (9) 

d t B, + —8 lu {m(V v ,B z -V t B w ))=0. (10) 
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2.2. Reduction to time-dependent ordinary differential 
equations 



Let us define the Jeans frequency: 



and 

t = t/SIj, 



(12) 



where po is the gas density. We consider the following 
fields: 



p(r, 


zu, 
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= Pod{r), 


P(r, 


ZU 
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= Po&j (P c (t) - Ptu(t)zu 2 - 


*(r, 


ZU, 
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= -r! 2 7 d(r)(A ro (r)w 2 +A z (t 


V„(r, 


ZU, 
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ZD 
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= 9( T )nj w, 
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z) 


- W) J ' 




ZU 


z) 


= V^oPo J h„(r) zu, 


Be(r, 


ZU, 


z) 


= VMoPo & J he(r) zu, 


B z (t, 


ZU 


z) 


= B c (t) + y/^opo fij h z (r) z. 



(13) 



where the dotted values denotes derivative against r and 
where d, P c , p m , p z , A m , A z , a, c, 6 , h^, h z , h g , and B c 
are functions of r. 

The definitions stated by Eqs. lfT3|l lead with Eqs. 
(|10|) to the following ordinary differential equations: 



Continuity: 

d(r) + (2^ + ^W)=0, 
V ayj) C ( T )J 

Energy conservation: 



(15) 
(16) 
(17) 



d(r)^l=2p z (r)-2A z (r)d(r) 2 , 
c(r) 



Nullity of magnetic divergence: 
(11) 2hUr) + h z (T)=0, 

Induction equation: 



Mr) - + 44^( r ) = °. 

a(r) c(r) 



3 

(21) 
(22) 
(23) 



2 44 + 44 ) W- *(r)(2/, x (r)+/,..(r)) = I I. ( 2-1- i 
a(r) c(r) 



a(r) c(r) 

y3 c (r) + 244 5 C M = 0- 
a(r) 



(25) 
(26) 



Eqs. lO, JEJl, (O, (113, 121, 123 and |23 are 

straightforwardly integrated and lead to: 



d(r) = 



1 



a(r) 2 c(r) ' 
o 



Pz(t) =P z -r-T2, 
c(t) 2 

hi 



(14) h e {r) -- 



-(P (r)d-^(r)) = 0, 

7(fc(r)<r(r)) + 2^(p ro (r)^(r)) = 0, 
ar a(r) 

A(p 2 ( T )rf-7 (r )) + 2^\(p z (T)d-^T)) = 0, 

dr c(r) 

Poisson equation: 
4A a (T)+2A z (r) = l, 

Momentum conservation: 

d ( T ) f44+^( T )) =2p~(T)-2A m (T)d(T) 2 ~2ho( T y,(19) 



h z (r) : 
Bc(r) 



a(r) 2 c(r) ' 
o,(t) 2 c(t) 



a(r) 2 c(r) 
5° 



a(r) 2 ' 



(27) 
(28) 
(29) 

(30) 

(31) 
(32) 
(33) 
(34) 



where P°, p%, p° z , h%, hg, B° are real numbers. 

Defining W(r) = a(r) 2 0(r), the system of Eqs. (|ri|) - 
(|26|l reduces after some algebra, to the three differential 
ordinary equations: 



o(r) 



(18) 5(t) = 



c(r) 



2p°, , l^ 2 (r) 2A ro (r) 2(^) 2 



a(r) 2 f- 1 c(r)T- 1 a(r) 3 a(r)c(r) o(t)c(t) 
2p° 2A z (r) 



,(35) 



Mr). 



d(r) ( 0(r) + 2^y^(r) ) = 2^(t)/i s (t), 



(20) W ^ = ^r) 



a(T) 2 -r- 2 c(Tp a (r) 2 ' 



(36) 
(37) 
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2.3. Gravitational potential 



2.4. Physical interpretations 



In the system of Eqs. (135(1 - (137(1 . the gravitational force (i.e. 
A ro and A z ) is unknown and must be specified. Its value 
indeed depends on the boundary conditions. For simplic- 
ity, it can be assumed that A^, and A z remain constant 
along time, as it was the case in the studies carried out 
by Hennebelle (2001). However the problem of boundary 
conditions has not been addressed for these solutions. It is 
likely that because of boundary conditions, rotation and 
magnetic field will induce an evolution of the shape of the 
cloud (and consequently a variation of A^ and A z ) for 
most of these solutions (if cylindrical geometry is assumed 
consistent solutions can be obtained, see Hennebelle 2003). 

More consistently, if one considers elliptical bound- 
ary conditions, i.e. a uniform axisymmetric spheroid, then 
one knows the gravitational potential as a function of 
the spheroid eccentricity e (Mestel 1963, Lin et al. 1965, 
Chandrasekhar 1969, Binney & Tremaine 1987). This al- 
lows to address properly the problem of boundary condi- 
tions and to obtain self-consistent solutions. 

For an oblate spheroid, we have: 

1 \/l — e 2 / arcsin(e) 
4^ 



Me) 



ivr 

2 e 



1 



- y/l -e 2 
arcsin(e) 



s/T 



(38) 







1 



and consequently: 

' A m - 1/6 
A z -> 1/6 

A m ~ (c/a) 2 
A z - 1/2. 

For a prolate spheroid, we have: 







(39) 



Me) 



1 1 
1~ 
1 1 
2 _ 



1 



1 -e 2 
1 



2e 



-In 



1 

~ 2e 
1 + e 
1 - e 



-hi 



1 - e 



(40) 



and: 







A v 
A z 

A z 
A, 



,1/6 
■1/6 

.1/4 

(a/c) 2 ln(a/c) 

lit, and l z being the length of the radial and the polar 
axis respectively, the eccentricity is given by: 







(41) 



e(r) = Wl 



for an oblate spheroid (l w > l z ), and by: 



e(r) 



(42) 



(43) 



for a prolate spheroid (l w <l z ). 

The lengths of the axis evolve with time according to: 

l„( T ) = k(0)a(r) , I z (t) = i,(0)c(r). (44) 



The system of Eqs. lJ35 )l -l|37 jl is of second order in a and c 
and of first order in W. It thus depends on 5 initial con- 
ditions, however one can assume without any restriction 
that a(0) = 1 and c(0) = 1. We thus have the 3 initial con- 
ditions, a(0), c(0), and W(0) which respectively represent 
the initial radial velocity, axial velocity and rotation. 

We also have the 5 parameters P°, 7, h%, hg and e(0) 
(as we will see with Eq. (JJ^J, e(0), P°, p% and p° z are in- 
deed related) which represent respectively the initial ther- 
mal pressure, the polytropic index, the initial magnetic 
poloidal field and magnetic toroidal field and the initial 
eccentricity. 

The system of Eqs. (13511 . (136(1 and 1(37(1 is much sim- 
pler than the system of Eqs. Q-lJTDJ and in spite of im- 
portant physical restrictions, allows to study various non- 
linear aspects of the gravitational contraction. Before solv- 
ing Eqs. (|3*5|l - (|3*7jl numerically, and in order to anticipate 
the numerical results, we discuss the physical meaning and 
mathematical behaviour of each of the terms in these three 
equations. 

2.4.1. The radial momentum conservation 

Eq. ((35(1 derives from the radial momentum equation, the 
first term of the right-hand side is the thermal pressure, 
the second is the centrifugal force, the third is the gravita- 
tional force, and the last is the toroidal magnetic pressure. 
The poloidal magnetic pressure does not appear since it 
cancels out (the poloidal magnetic field is force free) as 
it was already the case in the studies of Aburihan et al. 
(2001) and Hennebelle (2001). 

Since the density is uniform, the thermal pressure gra- 
dient is due to a non uniform temperature field. Although 
such temperature gradients are expected if the gas is 
opaque to its radiation, in the case of dense cores that we 
will consider in the next sections the temperature should 
be nearly constant. Rigorously speaking, these temper- 
ature gradients should be seen as ersatz that allows to 
treat the thermal pressure of a uniform density cloud semi- 
analytically. 

• If the collapse is spherical (p% = p z , e(0) = 0, 
W(0) = 0, hg = 0), a = c and Eq. can be integrated 
once leading in the case 7 ^ 1 to: 



1 , 2 
-a = 

2 3- 



Pi 



3 a 



(45) 



3 7 a 3 ^" 3 

where K is a constant which depends on the initial con- 
ditions. If 7 < 4/3 the collapse is possible, since the grav- 
itational force dominates the pressure force in the limit 
a — > 0. The collapse is not possible if 7 > 4/3. 

• If the cloud contracts radially, we can assume 
that c remains constant in the final stage. Remembering 
Eq. 1(41(1 . Eq. 1(35(1 can be further integrated (provided 
7 7^ 1) in the limit a <C cq (c(t) ~ Co): 

> p% 2M/ 2 (0) 



1 



27 a 2 t- 



■2„7" 



(46) 
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1 

2cq" 



Co 



ln(a) + K, 



where K again is an integration constant. Therefore, ther- 
mal pressure cannot support the cloud if 7 < 1 whereas if 
7 > 1, it prevents the collapse. The centrifugal force pre- 
vents the collapse and the Lorentz force due to the toroidal 
magnetic field adds up to gravity. 

• If the collapse occurs along the polar axis only 
(c — > 0), the pressure force is proportional to: c 1-7 , grav- 
ity is proportional to c (remembering Eqs. I|39|) . oc c 2 ) 
and the toroidal force is proportional to 1 / c. Consequently, 
the radial component of the gravitational force becomes 
small compared to the thermal pressure whereas on the 
contrary, the toroidal magnetic force dominates the ther- 
mal pressure. 

2.4.2. The axial momentum conservation 

Eq. H36J) derives from the axial momentum equation. It 
contains only the thermal pressure (first term of right- 
hand side) and gravity (second term) since the magnetic 
forces in our geometry are radial and azimuthal. 

• If the cloud contracts radially (a — > 0), then with 
Eqs. H41JI . we find that the pressure force is proportional 
to a 2 ~ 27 and gravity to ln(et), consequently the thermal 
pressure dominates if 7 > 1. 

• If the collapse occurs along the pole (c — ► 0), assum- 
ing that a remains constant (a ~ ao), we can integrate 
Eq. (|36|1 once (assuming that 7 ^ 1): 



1 



7 ag 7 - V- 1 



jC + K, 



(47) 



consequently, the collapse into an infinitely thin disk is 
not possible if 7 > 1. As pointed out in the conclusion of 
Lin et al. (1965), at some point the cloud may rebound. 

2.4.3. The azimuthal momentum conservation 

Eq. I|37ll derives from the azimuthal momentum equation 
(Eq. 0J|. The term of the right-hand side is the magnetic 
tension. If h^hg < 0, it decreases the rotation velocity 
(magnetic braking) and tends to increase it in the other 
case. 



2.4.4. Boundary conditions 

In this section, we address the problem of the boundary 
conditions. 



• Thermal pressure 
The cloud boundary {wb, Zb) is given by the condition that 
the thermal pressure vanishes. Eqs. (|13|1 lead to: 

P c (t) - v^l ~ vA = 0. (48) 

Thus, at the points (a(r)^ ro (0), 0) and (0, c(t)I z (0)), we 
have with Eqs. f5%|l - ll3U|) : 

F° = p% (^) 2 +P z (^) 2 =plU0f=p%(0f. (49) 



Therefore the boundary is a spheroid having axis lengths 
equal to: Z ro (0) = y/P^Jp% an l z (0) = ^P°/p° z . We have 
consequently: e(0) = \J\ — P^/p® for an oblate cloud and 
e(0) = \J\ — p z /p% for a prolate cloud. In the case of a 
sphere, one has: 1-^(0) = l z (0) implying p% = p° z . 

• Magnetic field 
The component of the magnetic field which is perpendic- 
ular to the cloud surface must be continuous through the 
surface whereas the component which is tangent to the 
cloud surface can jump discontinuously from a value to 
another leading to surface electric currents. 

Since the density of the external medium is zero, the 
Lorentz force must vanish in the external medium and the 
Alfven speed is infinite. 

This implies that the solutions given by Eqs. (|35[1 - 
(|37[1 with the condition stated by Eq. (|49[) are compat- 
ible with a poloidal magnetic field having the same struc- 
ture (Eqs. I13JI in the external medium and in the cloud 
(h^ ^ implies that the magnetic field becomes infinite 
when r = \J vo 2 + z 2 — * 00) and a vanishing toroidal mag- 
netic component in the external medium. Note that the 
structure of the magnetic field does not necessarily follow 
the self-similar forms stated by Eqs. (|13|) in the external 
medium and that other structures may be solution of this 
problem. 

The structure of the poloidal magnetic field in the ex- 
ternal medium adjusts instantaneously as the boundary of 
the cloud moves, so that the field in the external medium 
remains force-free, since the Alfven speed is infinite. 

2.4.5. Energy balance 

In the hydrodynamical case (h° e = h% = 0), the solution 
is defined in a finite domain of space, and it is possible to 
find the energy of the corresponding cloud. 

Let us multiply Eq. by 2Z ro (0) 2 d, Eq. GHl by 

l z (0) 2 c and add the two expressions. Integrating the re- 
sulting equation, we obtain when 7 7^ 1: 

1 (2Mo v + imv) + 7 ! la27 p 2 ° 7 _, + 



For 7=1, the integrated equation reads: 

i (2U0) 2 d 2 + l,(0)V) - 2P° ln(a J o) + 

umvH») ( iUor ^ +2Uof ^) =Ki (51) 



which is a first integral of Eqs. H35fl - i|3(j|) and represents 
the total energy of the system. The first term is the kinetic 
energy, the second one is the thermal energy, the third one 
is the rotation energy and the last one is the gravitational 
energy. 

If the magnetic field is non zero, the energy of the 
system is infinite since the magnetic energy of the sur- 
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rounding gas is infinite and it is not possible to find such 
a relation. 

The a and (3 parameters, respectively the ratio of the 
thermal to gravitational energy and of the rotational to 
gravitational energy are equal to: 



a(r) 



2P° 



,4-2 7 „2- 7 



1 -\{Al zo {0fA^ + 2l z {0) 2 A z c^y 

l m {0fW{0fc 
(4^(0) 2 A ro a 2 + 2/40) 2 A z c 2 )' 



(52) 



(53) 



In the case of an initially spherical cloud, we simply have: 
a(0) = 2P°/( 7 - I) and /3(0) = W(0) 2 . 

2.5. Potential and weaknesses of the approach 

Although the derivation of Eqs. JTJ-fTOJl into Eqs. 
l|37|) is exact, the resulting reduction leaves aside some 
important physical effects. Nevertheless, it allows to catch 
some complex non-linear behaviours of the collapse of a 
rotating and magnetized polytropic cloud at very little 
cost. This approach complements a full numerical simu- 
lation for which all the parameter space cannot be easily 
explored and which can possibly be affected by the nu- 
merical resolution and by the numerical scheme. Finally 
such simple solutions can be used to derive semi-analytic 
criteria (see e.g. Tsuribe & Inutsuka 1999). 

We now discuss the potential and the weaknesses of 
this reduction. 

2.5.1. Potential 

Since the model is two dimensional, it is well suited to 
the study of the global behaviour of a collapsing core in- 
cluding the radial and axial rebounds. It allows to study 
the effects of rotation on the collapse, including the shape 
evolution and the subsequent variation of the gravitational 
potential. 

Since the equation of state is polytropic, the model 
also allows to study the influence of the polytropic index 
on the cloud dynamics. 

Since it includes magnetic field, it allows to catch some 
magnetic effects as well, like toroidal pinching and mag- 
netic braking. 

2.5.2. Weaknesses 

The homologous reduction stated by Eqs. I|13[) excludes 
any complex hydrodynamical features like shocks and tur- 
bulence. Due to homology, the dissipation by viscosity 
vanishes. 

The temperature decreases from the center to the 
cloud edge. This prevents the inwards propagation of the 
rarefaction wave that is usually obtained in the numerical 
simulations of cloud collapse. 

The reduction stated by Eqs. I|13|) leads to a vanishing 
magnetic poloidal support. This clearly means that the 



solutions are not able to describe the cloud evolution if the 
poloidal magnetic support is the dominant process as it is 
the case for the collapse of a subcritical cloud controlled 
by ambipolar diffusion (e.g. Basu & Mouschovias 1994). 

In these solutions the growth (or the decrease) of the 
toroidal magnetic component due to the stretching of the 
poloidal component by differential (radial or axial) rota- 
tion vanishes, since the corresponding term cancels out 
in Eq. I|24|l. This effect however is potentially important 
since the collapse could be induced by the toroidal mag- 
netic pressure (Habe et al. 1991, Tomisaka 1991). The ax- 
ial support due to the toroidal magnetic field vanishes as 
well. 

3. Isothermal cloud collapse: phenomenology 

In this section we investigate the dynamical behaviour of 
the solutions of Eqs. (jTJ-JTni) described by Eqs. 
semi-analytically. The collapse of a cold cloud (with no 
thermal pressure) has been studied by Lin et al. (1965). 
They find that the oblate clouds collapse into a disk 
whereas the prolate clouds collapse into a spindle. This 
is a purely gravitational effect. The gravitational force is 
higher along the minor axis than along the major one and 
the initial asymmetry is amplified. 

3.1. Application to dense cores 

As an application we will consider the collapse of an 
isothermal prestellar cloud. The density profiles of prestel- 
lar clouds have been found to be rather flat in the central 
part (e.g. Bacmann et al. 2000, Motte & Andre 2001) hav- 
ing rj — dln[p]/dln[r] < 1. Consequently the uniform den- 
sity approximation should be acceptable to describe the 
dynamical evolution of their inner part during the precol- 
lapse phase. Once the collapse phase begins, strong density 
gradients develop, the present approach is no more valid 
and other analytical techniques must be used to describe 
the collapse phase (see e.g. the references quoted in the 
introduction) . 

The typical physical parameters of the prestellar cloud 
are: temperature ~ 10 K, density ~ 10 4 cm -3 (p ~ 3 10 -20 
g), size Iq ~ 0.1 pc, adiabatic index 7 = 1, rotation rang- 
ing from P = to j3 = 0.1, where (3 is the ratio of the 
rotational energy to the gravitational energy. The mag- 
netic field unit is y/nopo^ljlo — 25 pG. The time unit, 
is: ~ 2 10 5 years. 

3.1.1. Initial conditions 

As most of the observed cores are not spherical but have 
an observed aspect ratio of ~ 0.55, we will consider ini- 
tially oblate and prolate spheroids, starting with an initial 
aspect ratio around 0.5. 

For simplicity, we start with a thermal pressure such 
that the pressure force along the minor axis (along which 
it is stronger than along the major axis) is 3/4 of the 
gravitational force along the same axis. The pressure along 
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the major axis is then deduced from Eq. I|49l) . We start 
with initial velocity equal to zero, d(0) = and c(0) = 0. 

3.1.2. Numerical simulations 

In order to assess the validity of the solutions and under- 
stand in which regime they can be used, we perform hy- 
drodynamical numerical simulations using the SPH code 
described in Hennebelle et al. (2003). We use 50,000 parti- 
cles to describe the cloud and a constant external pressure 
equal to the initial pressure of the cloud. In the following, 
we simulate the evolution of clouds of initially uniform 
density. Note that this configuration is not identical to 
the configuration described by Eqs. 1)35(1 - 1137(1 . In the sim- 
ulation the rarefaction wave rapidly leads to non uniform 
density. Moreover the boundary pressure is not zero. 

In order to have a realistic comparison, the mass of the 
cloud is chosen equal to one solar mass and the tempera- 
ture to 10 K. 

3.2. Oblate cloud 

We first consider non-rotating (W(0) = 0), unmagnetized 
(q(Q) = kg I a(0) 2 / 'c(O) = 0), oblate clouds and then clouds 
with various rotation speeds. 

3.2.1. Non-rotating cloud 

We consider three different values of the aspect ratio, 
namely r(0) = Z*(0)/Z ro (0) = 0.7, 0.5 and 0.3. The nu- 
merical integration is stopped when the density reaches 
the value 10 5 . The results are displayed in Fig.^ the evo- 
lution of a(r), the radial axis (dotted line) and c(r) the 
polar axis (full line) are presented. 

• In the first case (r(Q) = 0.7), it is seen that, as in 
the study of Lin et al. (1965) gravity amplifies the initial 
anisotropy and the cloud becomes more and more oblate. 
It eventually collapses into a disk at t ~ 2.7/f2j. 

• In the second case, the minor axis collapses more 
rapidly than the major axis but after a strong contrac- 
tion (c(0)/c m in — 10 3 ), a rebound due to thermal pressure 
occurs at t ~ 2.3/Qj, the axial velocity becomes positive 
(outward motion) whereas the cloud still contracts radi- 
ally. Let us recall that in their conclusion, Lin et al. (1965) 
predict that the cloud should oscillate if thermal pressure 
would be introduced. The cloud finally collapses into a 
spindle at t = 2.6/f2j. Shortly before, the collapse along 
the polar axis restarts. Note that the rebound is very stiff 
and some care is required in the numerical resolution. We 
use adaptive time steps imposing that the maximum vari- 
ation of a and c between two steps is less than 1 % and 
we checked that the same result is obtained by requiring 
a variation of 0.5 %. We also check that the energy given 
by Eq. 1(51(1 is well conserved. 

• In the third case (r(0) = 0.3), the rebound oc- 
curs still at t ~ 2.6/f2j but the contraction is weaker 
(c(0)/c m j n — 1/30). This is due to the fact that since the 
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Fig. 1. Semi-analytical solutions of the collapse of an 
oblate isothermal (7 = 1) cloud with initial axis ratio 
equal to 0.7 (upper panel), 0.5 (middle panel) and 0.3 
(lower panel). The evolution of the cloud radial axis con- 
traction factor, a(r) (dotted line), and the cloud polar axis 
contraction factor, c(r) (full line), are displayed. 



cloud is initially more oblate than in the previous case, 
the axial velocity is smaller leading to a weaker contrac- 
tion. Also the dynamical time scale along the minor axis 
is shorter and the radial contraction is weaker when the 
rebound occurs and the density is smaller by the time of 
this rebound. After this rebound, the collapse along the 
polar axis restarts and finally, the cloud collapses into a 
disk, i.e. the density reaches our threshold, at t ~ 2.9/ttj. 

• Fig. [3 presents the numerical simulation of a non ro- 
tating oblate cloud (the initial aspect ratio equal is 0.5). 
The initial length of the major axis is 0.07 pc. Five time 
steps are displayed. It is seen that the velocity field (both 
radial and axial components) can be well approximated by 
an homologous law in the inner part and that the inner 
part is also reasonably flat during the precollapse phase. 
Because of the rarefaction wave, not describe by the semi- 
analytical solutions, the density contrast in the outer part 
increases and the velocity is no more homologous. It is 
seen that at time t = 2.23/f2j the cloud bounces axi- 
ally because of the thermal pressure. Note that the inner 
part slows down axially and stops collapsing (t = 1.78 
and 2.23/ ft j) without bouncing. The evolution of the ra- 
dial component, V m is monotonous in the inner and in 
the outer part. At time 2.67/ilj it is seen that the col- 
lapse along the polar axis has restarted in the inner part 
whereas the outer part is still bouncing. Although this 
behaviour is too complex to be completely described by 
the semi-analytical solutions, they nevertheless undoubt- 
edly reproduce some complex features of the cloud evolu- 
tion including the axial rebounds and subsequent collapse. 
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Fig. 2. Numerical simulation of a rotating oblate and 
isothermal (7 = 1) cloud with initial axis ratio equal 
to 0.5. First panel displays the equatorial density, second 
panel the density along the polar axis, third is the radial 
velocity at the equator and fourth is the axial velocity 
along the polar axis. Full line corresponds to t=0.89, dot- 
ted line to 1.34, dashed to 1.78, dot-dashed to 2.23 and 
triple-dot-dashed to 2.67 



Note that the largest density reached before the rebound 
is significantly overestimated by the semi-analytical so- 
lution. This is due to the thermal pressure which is not 
accurately calculated. However, the duration of this phase 
of strong condensation is short and does not alter much 
the subsequent evolution. 



Fig. 3. Semi-analytical solutions of the collapse of a rotat- 
ing oblate and isothermal (7=1) cloud with initial axis 
ratio equal to 0.5. Upper panel is W(0) = 0.2, middle one, 
W(0) = 0.3 and lower panel is W(0) = 0.4. The evolution 
of a(r) (dotted line) and c(r) (full line) are displayed. 



3.2.2. Rotating oblate cloud 

We consider the evolution of an initially oblate cloud with 
initial axis ratio r(0) = 0.5 and investigate three initial 
rotation values, W(0) = 0.2,0.3 and 0.4. The results are 
displayed in Fig. |31 

• The first case (VF(0) = 0.2) is similar to the evolution 
of the non-rotating clouds (see Fig. except that due 
to the centrifugal support, the rebound occurs later (t ~ 
3/f2j). The cloud collapses into a disk ati~3.3/fij. 

• In the second case (W^(0) = 0.3), it is found that the 
cloud does not collapse but makes a series of rebounds. 
This is due to the fact that rotation strongly supports 
the cloud and consequently the contraction along the po- 
lar axis is too weak (density is lower and gravitational 
force along the polar axis smaller than for the previous 
case) for allowing the density to reach our threshold. The 
cloud presents a quasi-periodic behaviour. Since Eqs. (|35|l - 
()36fl are non-dissipative, the cloud never relax and never 
collapse. In a more realistic situation, the number of re- 
bounds depend on the typical time scale of the dissipative 
processes. If this time is short compared to the oscillation 
period, the rebounds will be damped rapidly (see Fig. El 
for a numerical estimate of this time). 

• In the third case (1^(0) = 0.4), the centrifugal force 
is stronger than for the previous case and the oscillations 
have a longer period and a smaller amplitude. 

• Fig. ^presents the numerical simulation of a rotating 
(W(0) = 0.3) initially oblate core (r(0) = 0.5). The initial 
length of the major axis is 0.07 pc. As in the case studied 
in the previous section, it is seen that the homologous ap- 
proximation is fair to describe the velocity field and that 
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Fig. 4. Same as Fig. [3 for a rotating oblate cloud (initial 
aspect ratio equal to 0.5 and W(0) = 0.3). Full line cor- 
responds to t=0.89, dotted line to 1.78, dashed to 2.68, 
dot-dashed to 3.57 and triple-dot-dashed to t = 4.46/Qj. 



the cloud is rather flat in the inner part. The cloud starts 
to contract axially and radially. At time t = 2.68/flj the 
cloud rebounds axially and at time 3.57/S7j, the inner part 
rebounds radially whereas the outer part is still contract- 
ing. At time t = 4.46/Qj the outer part of the cloud ex- 
pands radially whereas the inner part is contracting. This 
behaviour is again very similar (though more complex) to 
the case displayed in the second panel of Fig. |3| where the 
axial rebound occurs at time ~ 3. 5/17 j and the radial one 
at ~ 4.2/Slj. The disagreement on the times is about 20 
%. 

During the subsequent evolution, the cloud continue 
to bounce axially and radially several times and looses its 
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Fig. 5. The kinetical energy as a function of time. Full 
line is the numerical solution and dotted line is the semi- 
analytical one. 



energy slowly as can be seen in Fig. [3] where the evolution 
of the kinematical energy of the cloud is displayed. The 
minima correspond to the most expanded states and the 
maxima are close to the most condensed. It is seen that 
the global behaviour of the cloud is reproduced by the 
semi- analytical solution during few oscillation cycles. The 
details however are not accurately described. In particular, 
as already mentioned, the semi-analytical solutions over- 
estimate the largest density reached during the oscillation 
cycles. 

It is worth noting that such an oscillating behaviour 
may have been recently observed in the dark cloud 
Barnard 68 by Lada et al. (2003) where both systematic 
inwards and outward motions have been inferred. 

3.3. Prolate cloud 

We consider now the collapse of prolate clouds, i.e. clouds 
having filamentary shape with an initial aspect ratio 
r(0) = i w (0)/Z z (0) = 0.5. We first investigate the case 
of rotating and unmagnetized clouds and then rotating 
and magnetized filaments. 

3.3.1. Rotating unmagnetized cloud 

Three values of the rotation speed are considered, W(Q) = 
0, 0.15 and 0.3. The results are displayed in Fig.|SJ 

• The first case (upper panel) is similar to the result of 
Lin et al. (1965). The major axis collapses first since grav- 
ity is stronger in this direction than along the major axis 
(A z /A m ~ 0.42) and the cloud collapses into a spindle. 

• In the second case (W(0) = 0.15), it is seen that 
the collapse of the minor axis is delayed by the centrifugal 
support and that the dynamics of the two axis is compa- 
rable. 

• When the cloud rotates more rapidly (W(0) = 0.3), 
the centrifugal force is stronger and the cloud collapses 
more slowly in the radial direction than along the polar 
axis. The cloud collapses into a disk at time t ~ 3.5/fij. 
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Fig. 6. Semi-analytical solutions of the collapse of a ro- 
tating prolate isothermal (7 = 1) cloud with initial axis 
ratio equal to r(0) = 0.5. Upper panel is W(0) = (no 
rotation), middle panel, W(Q) — 0.15 and lower panel, 
W(0) = 0.3. The evolution of a(r) (dotted line) and c(r) 
(full line) are displayed. 



Fig. 7. Semi-analytical solutions of the collapse of a ro- 
tating and magnetized prolate isothermal (7 = 1) cloud 
with initial axis ratio equal to r(0) = 0.5. Upper panel is 
W(0) = 0.1,9(0) = 0.3 middle panel, W(0) = 0.1, q(0) = 
0.7 and lower panel, W(0) = 0.3, q(0) = 0.7. The evolution 
of a(r) (dotted line) and c(r) (full line) are displayed. 



3.3.2. Rotating and magnetized cloud 

We investigate now the collapse of a rotating and mag- 
netized cloud. For simplicity, we consider the same initial 
condition as in Sect. 13.3.11 and add a toroidal magnetic 
field with an amplitude q = 0.3 and 0.7. The results are 
displayed in Fig. 

• In the first case (W(0) = 0.1, q(0) = 0.3), the 
evolution is similar to the non rotating prolate cloud 
(PF(0) = 0, q(0) — 0) but the minor axis collapses more 
rapidly. This is due to the toroidal magnetic field that 
adds up to gravity in the equatorial plane but not along 
the polar axis (initially we have, A z /(A m +q{0) 2 ) ~ 0.29). 
At time t = 1.9/Clj a radial rebound occurs because ro- 
tation prevents the collapse. The cloud then becomes less 
and less prolate, until t ~ 2.5/flj, where the collapse in 
the radial direction restarts. It eventually collapses into a 
disk at time t ~ 2.7/Clj. 

• If the toroidal magnetic field is stronger (q(0) = 0.7), 
the cloud is more strongly squeezed [A Z /(A^ + q(0) 2 ) ~ 
0.12) and the minor axis collapses more rapidly than in 
the previous case and at t ~ l.2/Qj, the first rebound 
occurs. By this time, the major axis has little evolved be- 
cause its dynamical time scale is slower. At t ~ 1.9/f2j, 
the collapse in the radial direction restarts and the cloud 
makes a second rebound at time t ~ 2.4/Oj. It eventually 
collapses into a disk at t ~ 2.5/fij. 

• The evolution of the more rapidly rotating cloud 
(W(0) = 0.3) is very similar to the slowly rotating one. 
The only difference appears when rotation becomes sig- 
nificant, i.e. during the contractions and the subsequent 



rebounds at t ~ 1.2 and 2.5/f2j. This is due to the fact 
that in these two cases the centrifugal force is small com- 
pared to the toroidal pinching when the cloud is weakly 
condensed. The centrifugal force plays an important role 
only when the cloud is significantly radially condensed. 

Let us summarise our conclusions. The collapse of a 
prolate and magnetized core (with strong toroidal pinch- 
ing) occurs slowly through gravity along the polar axis and 
rapidly through toroidal pinching and gravity in the ra- 
dial direction. Since rotation prevents radial collapse, few 
rebounds occur during the slow collapse along the polar 
axis and finally the cloud collapses into a disk. Because 
of the toroidal magnetic field, the dynamical time scale of 
the minor axis is much shorter than the dynamical time 
scale of the major one. This process depends weakly on 
the rotation amplitude (as long as it is neither zero nor 
very large) and on the thermal pressure. 

4. Apparent shape of the clouds, comparison with 
observation 

In spite of the important restrictions listed in sect. I2.5~21 
the solutions presented in this paper enable the study of 
the cloud shape evolution during the precollapse phase. It 
is worth trying to compare the derived predictions with 
observations of dense cores traced by the NH3 molecule 
(Jijina et al. 1999, Lee & Myers 1999). Since our simple 
model assumes a uniform density and does not include 
complex processes like accretion or protostellar outflows, 
we select the starless dense cores only from the Jijina et 
al. (1999) catalog. 
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Two recent studies (Jones et al. 2001, Goodwin et al. 
2001) analyse the observed apparent axis ratio distribu- 
tion derived from this catalog and conclude that this dis- 
tribution is compatible with dense cores being triaxal. 

Since our model makes restrictive assumptions and 
since the Jijina et al. (1999) catalog is not an homoge- 
neous sample, the cores having been observed with var- 
ious instruments and at various resolutions, the present 
analysis is rather tentative. The idea is to see if a good fit 
can be obtained and for which range of parameters. The 
result will not be more than indicative. 

4.1. Method 

The first step of our analysis is to derive a shape distribu- 
tion from the time sequence of a and c. This is achieved by 
counting the fraction of time that one cloud spends with 
a given axis ratio. 

The second step is to project the axis ratio distribution 
along the line of sight. Let ip be the axis ratio distribution, 
and 0, the apparent (or observed) axis ratio distribution. 
The projection along the line of sight is achieved according 
to the law (Binney 1978, Fall & Freck 1983): 

<P(r)=r [ ds(l- s 2 )- 1 ' 2 ^ 2 - s 2 )- 1 ' 2 ^) (54) 
Jo 

for an oblate distribution, 

4>{r) = r- 2 ^ dss 2 (l - s 2 )- l ' 2 {r 2 - s 2 )- 1 ' 2 ^) (55) 
Jo 

for a prolate distribution. 

The third step is to calculate the x 2 between this the- 
oretical distribution and the distribution derived from the 
catalog of Jinina et al. (1999) which includes 79 starless 
cores. 

We explore a large set of parameters: 

- the ratio between the 2 axis, r(0) = ^ 2 (0)/7 ro (0) varies 

from 0.5 (oblate cloud) to 2.2 (prolate cloud), with a 
regular sampling of 0.1 leading to 18 values, 

- the magnetic toroidal intensity q(0) varies from to 1.5 

(16 values), in physical units it ranges from to 40^iG. 

- the rotation velocity, 1^(0), from to 0.6 (7 values) 

which corresponds to (3 ranging from to 36%. 

- the thermal pressure from = (Ai(Q))/10 to Ai(Q) (10 

values) where i = z for oblate cores and i = w for 
prolate cores. 

Consequently, we solve 18 x 16 x 7 x 10 (~ 2 10 4 ) times 
Eqs. 135|l - 136|l and calculate the corresponding apparent 
axis ratio and the x 2 with the observed distribution for 
each of these solutions. 

4.2. Result 

The best agreement obtained is for Z OT (0) / Z z (0) = 0.5, q = 
0.7, W(0) = 0.3, P c = A ra (0)/7, which is very close to the 
second and third cases we have considered in Sect. 13.3.21 
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Fig. 8. Apparent axis ratio during the cloud evolution 
(upper panel), axis ratio distribution (middle panel) and 
apparent axis ratio distribution (lower panel). The dot- 
ted line (second panel) is the fraction of prolate core, 
the dashed line, the fraction of oblate cores. In the third 
panel the dot-dashed curve is the apparent axis distribu- 
tion of starless cores belonging to the catalog of Jijina et 
al. (1999) whereas the full line is the theoretical distribu- 
tion obtained for the initial conditions / ro (0)/^(0) = 0.5, 
q = 0.7 (magnetic toroidal intensity), W(0) = 0.3 (rota- 
tion) and p% = ^4 ro (0)/7 (thermal pressure). The value of 
X 2 is 0.56. 



(prolate magnetized cores). The theoretical (full line) and 
observational (dot-dashed line) arc displayed in the third 
panel of Fig. |H1 The agreement is good and better than 
in the study of Jones et al. (2001) and Goodwin et al. 
(2002). The corresponding value of x 2 is: 0.56 (assuming 
an observational Poisson noise and no theoretical noise). 
The distribution depends weakly on the value of rotation 
and thermal pressure and is much sensitive to the initial 
axis ratio r(0) = 1^(0) /l z (0) and to the magnetic toroidal 
intensity. 

We find 179 cases with x < 1 (87 have x < 0-8) similar 
to the best case. They have Z z (0)/^ ro (0) equal to 2 to 2.2, 
q equal to 0.6 to 0.9 and all possible values of W(0) and 

Pc 
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We also find 59 cases that belong to another family. 
They are oblate (i z (0)/i ro (0) = 0.5) fast rotating cores 
(W"(0) > 0.4) similar to the cases presented in Sect. I3.2T21 
However, the smallest x is 0.85 and most of them have 
X > 0.9. The most important disagreement is due to the 
fact that too many clouds with apparent axis ratio close 
to 1, are found. 

4.3. Discussion 

The apparent axis ratio distribution derived from the evo- 
lution of the prolate magnetized filaments agrees with the 
observed distribution. 

The rotating oblate cores agree less well with the obser- 
vational apparent axis ratio distribution. Indeed the core 
remains oblate and Ryden (1996) shows that an oblate 
core distribution cannot explained the observed apparent 
axis ratio distribution. 

We would like to emphasize the fact that the axis ra- 
tio distribution depends weakly on the initial rotation and 
thermal pressure values (almost all values of these param- 
eters are compatible with the observed axis ratio distri- 
bution) and strongly on the initial axis ratio and on the 
toroidal magnetic field. Consequently, the agreement is ob- 
tained by adjusting two parameters only whereas there are 
9 bins in the fitted distribution. Therefore the fact that we 
obtain a parameter range that fits the data is not a trivial 
consequence of a large number of parameters. Indeed there 
would not be acceptable fits without a toroidal magnetic 
field. 

The most striking feature of the unprojected axis ratio 
distribution of the example displayed in Fig. |H1 is that it 
is bimodal (with two peaks around p± ~ 0.5 and P2 — 
0.7). This is due to the fact that the cloud spends a long 
time, with its initial aspect ratio and a long time with the 
aspect ratio it has after the first rebound. This feature is 
smoothed by the projection effect and the apparent axis 
ratio distribution is only slightly or marginally bimodal. 

The data of Jijina et al. (1999) are compatible with 
(though marginally statistically significant) a starless 
cores apparent axis ratio distribution being bimodal (see 
the dot-dashed curve of Fig. |H1 or the starless cores axis 
ratio distribution showed in Goodwin et al. 2002). The 
same trend appears also in two of the subsamples studied 
by Ryden (1996) and in the data studied recently by Curry 
(2002) although in all of these samples the poor statistic 
precludes a definite answer. More data are needed to con- 
firm or to reject this feature. 

5. Conclusion 

New solutions of the gravo-magnetic contraction have 
been presented. They describe the evolution of a uniform 
density cloud with spheroidal shape and can be applied 
in various situations. These solutions allow to explore a 
wide set of parameters and constitute a useful comple- 
ment to numerical simulations. They can also be used as 
benchmarks in problems involving gravitational collapse. 



We have explored the phenomenology of the gravita- 
tional contraction of isothermal dense cores described by 
these semi-analytical solutions with special emphasis on 
the axial and radial rebounds. Numerical simulations have 
been performed to compare these solutions with the evo- 
lution of an isothermal initially uniform cloud and it has 
been found that they can reproduce some aspects of cloud 
evolution during the precollapse phase. The most impor- 
tant disagreement is the cloud outer part which is signif- 
icantly influenced by the rarefaction wave not taken into 
account in the semi-analytical solutions. 

For each of a large set of initial values, we have calcu- 
lated the apparent axis ratio distribution resulting from 
the evolution of the core and tentatively, compared it to 
the observational one derived from the catalog of Jijina 
et al. (1999). We find that the axis ratio distribution of 
a collapsing rotating, prolate (initial aspect ratio ~ 0.5) 
magnetized (~ 17 — 20/iG for a density of ~ 10 4 cm -3 ) fil- 
ament is in good agreement with the observational distri- 
bution (see Fig. [SJ . Definite conclusion cannot be drawn 
at this stage however since our model is very simplified 
and (mhd) numerical studies are required to confirm the 
results. Also the observational core sample is not homoge- 
neous and the statistical significance remains poor. Finally 
alternative distributions of evolving or static cores can fit 
the observational data as well. Nevertheless, we believe 
that prolate cores permeated by an helical magnetic field 
deserve further investigations. 
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